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Abstract 



We derive the discrete linear systems associated to multi-matrix models, the 
corresponding discrete hierarchies and the appropriate coupling conditions. We also 
obtain the Wi+oo constraints on the partition function. We then apply to multi- 
matrix models the technique, developed in previous papers, of extracting hierarchies 
of differential equations from lattice ones without passing through a continuum limit. 
In a q-matrix model we find 2q coupled differential systems. The corresponding 
differential hierarchies are particular versions of the KP hierarchy. We show that the 
multi-matrix partition function is a r-function of these hierarchies. We discuss a 
few examples in the dispersionless limit. 
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1 Introduction 



In 1-matrix models all the information is encoded in the Jacobi matrix Q, a semi- 
infinite matrix formed by three non-vanishing diagonal lines. This matrix remains the 
same in form whatever the polynomial potential is. This accounts for the relatively 
simple structure of 1-matrix models. When we turn our attention to multi-matrix 
models, the situation drastically changes. First of all we have several Jacobi matrices 
characterizing the same model. They are semi-infinite matrices with a finite or 
infinite band of non-vanishing diagonal lines. What is more important, the band size 
depends on the order of the polynomial potentials. 

As is well-known the purpose of all the matrix model analysis is to extract flow 
equations (w.r.t. to the perturbation parameters) for the partition function, as well as 
constraints known as string equations. It is believed that multi-matrix models should 
give rise in particular to generalized KdV-hierarchies |j^] which in turn relate to 
generalized Kontsevich models §, |],|i,|],§,@. 

As will be shown below, it is not difficult to extract lattice hierarchies and lattice 
string equations from multi-matrix models. In fact in a q-matrix model we find 
2q linear systems and, consequently, 2q lattice hierarchies. They are generalizations 
of the Toda lattice hierarchy and they are sometimes referred to as discrete KP 
hierarchies. However significant information can be extracted only from differential 
equations. The usual way to obtain a differential hierarchy from a lattice hierarchy is 
via a continuum limit. However, due to the aforementioned features of multi-matrix 
models, the application of the continuum limit technique has not proven as successful 
as in 1-matrix models, and the research in this subject has remained at a rather 
conjectural level. The only exception are the results obtained in concerning 
the VF-type constraints satisfied by multi-matrix models, and the analysis of some 
particular cases ||l^,||ll[| (see also H]). 

In a previous paper |]l^] (see also [0]) we introduced a method to extract a dif- 
ferential hierarchy from a lattice one. It consists essentially in using the first flow 
equation to eliminate the difference operations in the remaining equations. In this 
way the latter become differential equations with respect to the first flow parameter 
(which is interpreted as the 'space' coordinate). This procedure can be generalized 
to the lattice systems specific of multi-matrix models, and this is what we do in this 
paper. As our main result, we associate to any q-matrix model a set of differential 
hierarchies of flow equations, together with the appropriate string equations. The 
differential hierarchies turn out to be particular versions (reductions) of the general- 
ized KP hierarchy |15|. We will refer to these hierarchies as KP-type hierarchies or, 
when no confusion is possible, KP hierarchies for short. As we will see, the structure 
of multi-matrix models is much richer than it is usually believed. In particular it 
contains much more than the n-th KdV hierarchies. 

We recall that the importance of avoiding a continuum limit lies in the fact |14] 
that the differential hierarchies and string equations so obtained describe properties 
of the matrix-model lattice and therefore, presumably, topological properties. 

The matter is arranged as follows. Section 2 is a review of multi-matrix mod- 
els and an occasion to state our notations. In section 3 we represent an arbitrary 
multi-matrix model by means of suitable coupled discrete linear systems, and em- 
phasize the importance of the coupling conditions; the corresponding consistency 
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conditions result in (generalized) 2-dimensional Toda lattice hierarchies and string 
equations. In section 4 we work out in detail the Wi+oo constraints on the parti- 
tion function for two-matrix models. In section 5 we pass from lattice to differential 
formalism and find the KP structure of the differential hierarchies pertinent to multi- 
matrix models. We also prove that the multi-matrix model partition function is in 
fact the r function of the appropriate KP-type hierarchy. In section 5 we see a few 
explicit examples of the above differential hierarchy in the dispersionless limit and 
prepare the ground for future developments. 

2 Multi— Matrix Models: Partition Function 

The partition function of the q-matrix model is given by 
Z^it, c)= J dMidM2 . . . dMge^''^ 
where Mi , . . . , Mq are Hermitian N x N matrices and 

q q-l 
a=l a=l 

with potentials 

Pa 

K. = $]W^a a = l,2...,q. 

r=l 

The paS are finite or infinite positive integers. 

The ordinary procedure to calculate the partition function consists of three steps 

(i) . One integrates out the angular parts such that only the integrations over the 
eigenvalues are left, 

. q N 

Z^(t,c) = const j YlY[dXa,iA{Xi)e^A{Xq), 

a=l i=l 

where 

q N q-l N 

a=l i=l a=l i=l 

and A(Ai) and A(Ag) are Vandermonde determinants. 

(ii) . One introduces the orthogonal polynomials 

^n(Ai) = A" + lower powers, fln{\) = + lower powers 

which satisfy the orthogonality relations 

j dXi... dA,Cn(Ai)e^i(^i)+'^+^^(^')r/„(Ag) = K{t, c)6nm (2.1) 
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where 



q-l oo g-1 
a=2r=l Q!=l 

(iii). We expand the Vandermonde determinants in terms of these orthogonal polyno- 
mials and using the orthogonality relation ( p.lD , we can easily calculate the partition 
function 

N-l 

ZN{t,c) = const N\ Jl hi (2.2) 

1=0 

Like in 1-matrix models, knowing the partition function means knowing the coeffi- 
cients hn{t,c)^s. Like in 1-matrix models, we will show in the next section that the 
information concerning the latter can be encoded in suitable linear systems, which 
have the advantage of clearly showing the integrable character of the multi-matrix 
models. 



To end this section let us introduce some convenient notations. For any matrix 
M, we define 

{M)^j= Mij^, Mij = Mji, Mi{j) = M^j_i. 

We will say that the element Mi{j) belongs to the j-th sector. As usual we also 
introduce a natural gradation 

deg[Eij] =j-i. 

For a given matrix M, if all its non-zero elements have degrees in the interval [a, b], 
then we will simply write: M G [a, b]. 

Moreover M_|_ will denote the upper triangular part of M (including the main 
diagonal), while M_ = M — M_|_. We will write 

N-l ^ oo 

T¥(M) =Y.Mii, T>(M) = Ma 

1=0 1=0 

Analogously, if L is a pseudo differential operator, L+ means the purely differential 
part of it, while L_ = L — L^. 



3 The Coupled Discrete Linear Systems in q— 
Matrix Model 

In order to extract the discrete linear systems that characterize our q-matrix model 
we have to pass through some preliminaries. First we redefine the orthogonal poly- 
nomials in the following way 

^„(Ai) = e^^(^^)en(Ai), $„(A,) = e^''(^^)r?„(A,). 
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As usual we denote the semi-infinite column vectors with components ^'o, ^i, ■ ■ ■ , 
and ^2) • • • ) by ^' and respectively. With these polynomials, the orthogo- 

nality relation ( |2.lD becomes 

/ n d\p^n{\l)e^'^m{\q) = 6nmhn{t, c). (3.1) 
13=1 

This orthogonality relation is going to play a crucial role in our analysis. Next we 
introduce the following Q-type matrices 

Yl dXa"^ n{Xl)e^ K'^miK) = Qnm{a)hm = Qmn{a)hn, a=l,...,q. (3.2) 

a=l 

Among them, Q{l),Q{q) are Jacobi matrices: their pure upper triangular part is 
/+ = J2i Notice also that if we call 

Qn+l,niQ) = Rn+l 

we immediately find 

hn+l = Rn+lhn 

Therefore, like in 1-matrix models, we can express the partition functions in terms 
of Rn- 

Beside the above Q matrices, we will need two P-type matrices, defined by 

j n dK[-Q^-^n{\l))e^^m{\) = Pnm{l)hm (3.3) 



f d 

J dXi... dXg^n{Xl)e>' {-Q^^miXg)) = Pmn{q)hn (3.4) 



3.1 The coupling conditions or string equations. 



One important preliminary remark is that the matrices ( |3.2| ) we introduced above are 
not completely independent. More precisely all the (5(a)'s can be expressed in terms 
of only one of them and one matrix P. Expressing the trivial fact that the integral 
of the total derivative of the integrand in eq.(^) with respect to XaA < a < q 
vanishes, we can easily derive the constraints or coupling conditions 

P(l)+ci2Q(2) = 0, (3.5a) 
Co,-i,aQ{a - 1) + + Ca,a+iQ{a + 1) = 0, 2 < a < q - 1, (3.5b) 
Cq.i,qQ{q-l)+V{q)=Q. (3.5c) 

where we use the notation 

Pa 

Vi = Y.rta,rQ''-Ha), a = l,2,...,q 

r=l 

These conditions play an extremely important role in the study of multi-matrix 
models. We will also refer to them as string equations. 

Before proceeding further, let us pause a bit to make a few comments. 
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It is just these coupling conditions that lead to the famous Tyi+oo^constraints 
on the partition function at the discrete level [^](see section 4). 

These conditions explicitly show that the Jacobi matrices depend on the choice 
of the potentials. We can immediately see that these coupling conditions com- 
pletely determine the degrees of the matrices Q{a). Since P{1) — V{ and 
P{q) — Vq are purely lower triangular matrices, a simple calculation shows that 



where 



(5(a) G a =1,2, 



mi = {pq-l)... {p3 - l){p2 - 1) 

rria = (Pq - l)(p<?-i - 1) . . . {pa+1 -1), 2 < a < q - 1 

m„ = 1 



and 



TT-l = 1 

Ua = {Pa-1 - 1) ... (P2 - -1), 2 < a < q - I 

riq = {Pq~l - 1) ... (P2 - - 1) 

These equations show that, if we want a finite band structure for the Q{a), only 
a finite number of perturbations is allowed. Conversely, if we want, say, Q{1) 
to possess the full discrete KP structure - we should let pa — > oo for some a; 
then all these matrices will have infinitely many non-zero diagonal lines. So, 
generally speaking, these coupling conditions reveal the main difi^erence between 
multi-matrix models and one-matrix models. 

Finally these conditions suggest that there may be a relation between matrix 
models and topological models. As an illustration of this remark, we consider a 



particular case - the two-matrix models. In this situation, the eqs.(3.5a-3.5c) 
reduce to 

P(l) + ci2Q(2) =0, ci2Q(l)+P(2) =0, 

In particular 

ci2Q+(2) = -VliQil)), ci2Q+(l) = -V[{Q{2)) (3.6) 

Now, we have pi + P2 + 1 series of coordinates, i.e. Pi + P2 + 1 non-vanishing 
diagonals of Q{1) and Q{2). The above equations (|3.6| ) give {pi + P2) series of 
constraints. Therefore we are only left with one-series of independent coordi- 
nates, i.e. one can express all the elements of Q{1) and Q{2) in terms of the 
i?j's (and some coupling constants). This is not surprising since all the infor- 
mation of the matrix model is encoded in the partition function, which only 
involves Ri's. Finally the remaining coupling conditions provides more series of 
constraints, so that eventually there are no 'local' degrees of freedom left. This 
means that, if multi-matrix models admit any field theory interpretation, they 
should correspond to topological field theories. 
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3.2 The associated discrete linear systems 



The derivation of the hnear systems associated to our matrix model is very simple. 
We take the derivatives of eqs.( |3.1| ) with respect to the time parameters fa,r; and use 
eqs.(p.2|). We get in this way the time evolution of ^ and ^, which can be represented 
in two different ways: 

(*) Discrete Linear System I: 



d 



l<k<pi, 
l<k<pa; 



2<a<q, 



(3.7) 



[ |^M/(Ai) = P(l)^(Ai). 
The corresponding consistency conditions are 



1 



^ Q(i) = [Q(i 



d 
dt 



P(l) 



(3 
(3. 

(3.8c) 



( ^l8b|) and (|3l8^) are hierarchies, leferred to henceforth as discrete KP hierarchies, 
whose integrability and meaning will be discussed later on. The first equation ( f3.8a ) 
is often referred to as string equation. Throughout the paper we use the latter term 
also as synonym of coupling conditions (see subsection 3.1). 

(**) Discrete Linear System H: 



i>(A,) = Qiiqmxg) 



dt 



dt. 



$(A, 



Q^(a)$(Ag), l<k<p^; l<a<q-l 



(3.9) 



gJ-$(A,) = P(g)cI>(A,). 



with consistency conditions 



[Q{q), P{q)] = l, 

Ota,k 

_d 

dtah 



P{q) = [P{q), Q'Lia)] 



(3.10a) 
(3.10b) 

(3.10c) 



(***) The Other Discrete Linear Systems: 
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One may wonder what are the equations of motion of the matrices Q(a)'s, with 
a 7^ 1, g. In the Ught of the coupUng conditions it is not surprising to find out that 
we have enough information to write them down. They can actuahy be extracted 
from the eqs.( 3.8a - 3.8c ) (or ( |3.10a - 3.10c )) and eqs.( p.5a - 3.5c ). We have 



d 



d 



dt 



I3,k 



Q{a) = [Qi{P), Q{a)], 
Q{a) = [Q{a), Q^(/3)], 



1 < /3 < a 
a<(3<q 



(3.11a) 
(3.11b) 



where k runs from 1 to pp. 

It is tedious, but straightforward to prove that all the flows we have introduced 
so far commute with one another. 

It is natural at this point to ask whether eqs.( ^.lia . 3.11b| ) can be considered as 
integrabihty conditions of some hnear systems. The answer is yes, but the relevant 
linear systems are a bit different from the ones we met so far. Let us start from the 
orthogonal relation (|3.1|) once again, and define two series of functions 



(3.12) 



and 



4-)(t,A«)^ /'ndA^^„(Ai)e^-. 

/3=1 



(3.13) 



where we denote 



/3=a+l 



a— I oo a— I 

= X! X! */3,r"^/3 + C/3,/3+1 A/3 A/3+1. 
/3=2r=l /3=1 
q-1 oo g-1 

= X! X! ^I3,r^p + C/3,/3+1 A/3 A/3+1. 
/3=a+l r=l f3=a 



Obviously we have 

el^)(t,Ai) = ^n(Ai), 



ryW(t,A,) = cI>^(A,). 



but for other values of a one sees immediately that ^^"^ and are not polynomials 
anymore. However they still satisfy an orthogonality relation 



dA„el"Ht,A„)e^-(^")??i-)(t, A„) = 6nmhnit,c) 



VI < a < 



(3.14) 



The q series of functions £,^\t,Xa) span q spaces Ha- The q series of functions 
r/l"^ (t, Xa) also span q spaces 7^* 's, which are dual to Ha via the orthogonal relations 



t Tills structure bears a perhaps not accidental similarity with the structure of topological field 
theories IpI. 
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Now we extract the spectral equations and the time evolution of these new func- 
tions. From eqs.( |3.2| ), we immediately see that 



This tells us that the spectral equations are 

Aa^^") = Q(a)e("\ l<a<q. 
A„r/(") = Q(a)r/("), I < a < q. 



l<a< 



(3.15) 



(3.16) 
(3.17) 



On the other hand, from the definitions ( |3.12D and ( |3.13 ), and making use of eqs.( |3.7D 
and (^.9|), we can derive the evolution equations of the functions ^^^^ and r/("\ which 
take the following form 



d 



dtp^r 

d 



dt 



I3,r 



= -Q-_(/?)^(-), 



and 



d 



dtj3^r 

d 



dt 



-rj 



(a) 



f3,r 



(a) 



l<P<a-l, 
a<P<q. 

a + l<P<q, 
1 < /3 < Q. 



(3.18) 
(3.19) 

(3.20) 
(3.21) 



Now equations ( 3.16 - |3^2l| ) represent 2q discrete linear systems, two of which 
coincide with ( |3.7D and ( |3.9| ), the other {2q — 2) being the additional discrete linear 
systems which give rise to the compatibility conditions (|3.11a| , 3.11b| ) and to the 
appropriate string equations. 

Finally we notice that, thanks to the generalized systems of orthogonal functions 
just introduced, we can represent all the linear systems in a unified form as follows 



and 



Q{a)r]{Xa) = XalliXa), 
4^r?(A„) = Q;(/3)??(A„), « + !</?<, 
a 



(3.22) 



(3.23) 



dt. 



13,' 



??(Ac 



-QL{a)7]{Xa), l<P<a. 



One could phrase the above results by saying that introducing the functions 
ei"^(t,A„)'s and vt\t,X Q.) can be considered as a formal reductions of the multi- 
matrix models to 1-matrix models, since we are left with the orthogonal relations 
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(P^ , which are similar to the ones in 1-matrix models. The major difference comes 
from the fact that (,ii^\t, Aa)'s and ri>n\t, Aq) are very complicated functions rather 
than the polynomials 0. 

A comment is in order concerning the plethora of linear systems we found. First 
we notice that, due to the coupling conditions (3.7), any Q{a) can be expressed in 
terms of (3(1) and P(l) or, alternatively, of Q{q) and P{q)- This means that the 
linear system I or the linear system II can be considered as the fundamental ones. 
We already noticed that the equations of motion for any Q{a) can be derived from 
the flow equations of these two systems. In other words these equations of motion 
admit 2q different linearizations. In the following we will mostly be dealing with the 
systems / and II. However one cannot exclude that a linear system different from / 
and // might be more convenient for specific purposes. 

3.2.1 Coordinatization of the Q and P matrices 

Henceforth we will use the following coordinatization of the Jacobi matrices 

nil ™9 

Q{l)=U + J2J2''i(^^^,i-i^ Q(g) =/+ + ^^6Kx)Sm-« (3.24) 

i 1=0 i 1=0 

and for the supplementary matrices 



One should never forget that 



'3(«)=E E Ti''^^)Ei,i-U 2<a<q-l. (3.25) 

i l=~na 



ai{i) = I > i 

bi{i^=0 l>i (3.26) 



r/")(i)=0 l>i 



Let us start from the coordinatization (|3.24| ). We can see immediately that 

(^Q+(l)) 6j^i+i + ao{i)6ij, (^Q_{q)j^^= Ri6j^i-i 
Therefore we can write down the and tg^i-flows explicitly 

^^a/(j) = ai+i{j + 1) - az+i(j) + ai{j)(^ao{j) - ao{j - O) (3.27a) 

^ai{j) = Rj-i+iai_i{j) - Rjai_i{j - 1) (3.27b) 

Otq^l 

-T^hiij) = Rj^i+ik^iij) - Rjbi^iij - 1) (3.27c) 

-kU) = bi+i{j + 1) - 6;+i(j) + bi{j){bo{j) - boU - I)) (3.27d) 



dtg 



^ tI-\j) = Ti^lij + 1) - Tlllij) + Ti^\3)[a,{j) - a,{j - I)) (3.27e) 
Ti^\j) = R,^i+,tI%) - R,tI% - 1) (3.27f) 



"''Does this mean that multi-matrix models can be realized as 1-matrix models if we choose suitable 
non-polynomial interactions? 
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We will also need the following coordinatization of P{1) and P{q) 

ni) = E E (3.28) 

i oo 

Sq 



m=i: E Plm^,^-l (3.29) 



i l=—oo 

where si and Sg are suitable positive (possibly infinite) integers. 



3.3 Flow equations for the partition function. 

As a byproduct of the above results we can easily get the flow equations for the 
partition function 

lnZ7v(t,c) = Tr(Q'-(a)), (3.30) 



1 < r < Pa, I < a < q 

Using the consistency conditions, we can rewrite these equations in the following 
form 

lnZ^(t,c) = (Q^(a))^^_^, (3.31) 



1 < r < Pa, l<a< 



or equivalently 



lnZ^(t,c) = (Q'{a))^^ , (3.32) 

Otq^lOta,r ^ / N,N-1 

1 < r < Pa] I < a < q 

We end this subsection with one more remark. The information concerning multi- 
matrix models is encoded in the discrete hierarchy ( p.lla - 3.11bD . Once we solve 



these equations, we can reconstruct the partition functions of the matrix models: 
they will automatically satisfy ( 3. SO] ). 



3.4 The hidden 2-dimensional Toda lattice 

In subsection 3.2 we have shown that multi-matrix models are equivalent to cer- 
tain discrete linear systems subject to some constraints. If we ignore the coupling 
conditions, the discrete linear systems lead to discrete hierarchies. In fact we will 
show that these discrete or lattice hierarchies, in the q = 2 case, are nothing but 



2-dimensional Toda lattice hierarchies |21]||2^, and in the generic q case, they are 



2-dimensional Toda lattices with additional flows. 

Let us start from the first flow equations written above. In particular, from 



eq.( 3.27a ), we have 

d 



ao{j) = R,+i - Rj (3.33) 
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while, from( ^^ ), one can show that 

d 



dti,i 



Rj = RjiaoU) - ao{j - 1)) (3.34) 



Combining these two equations we obtain the following 2-dimensional Toda lattice 
equation 

-——InRj = Rj+i - 2Rj + Rj_i (3.35) 

Oti^lOtq^l 

In fact in terms of the coordinates (pj = Inhj, the above equation takes the familiar 
form 
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We could have gotten the same result from 

d 



Tj+l-Vj _ pVj-<Pj-l 



(3.36) 



and from 

d 



boU) = R,+i - Rj (3.37) 



dtq^l 

the latter being obtained again from (3.30). 



R, = Rj{bo{j) - bo{j - I)) (3.38) 



( 3.36| ) is nothing but the two-dimensional Toda lattice equation. In the q = 2 case 



(i.e. 2-matrix model) the other time flows generate the complete 2~dimensional Toda 
lattice hierarchy [^]. From the analysis in subsection 3.2, we know that the system 
is restricted by the coupling conditions, therefore we can say that 2-matrix models 
are nothing but constrained 2-dimensional Toda lattices. The q-matrix models with 
q > 2 contains not only the Toda lattice hierarchy, but also additional series of flows. 
It is perhaps interesting to remark that, from eqs. (|3.35 | ) and ( 3.36 ), the parameters 



ti^i and tq^i could be interpreted as space coordinates while the other parameters are 
really time flow parameters. Therefore the multi-matrix models hierarchies could 
underlie a 2+1 D theory. 

3.5 Integrability 

The discrete linear systems ensuing from multi-matrix models are integrable as is 
guaranteed by the fact that all the flows commute with one another. Most integrable 
systems we know possess a bi-hamiltonian structure with Hamiltonians in involu- 
tions. For the linear systems of multi-matrix models we are discussing here we can 
define a bi-hamiltonian structure, i.e. two compatible Poisson brackets. 

Let us be more explicit. Following [^], let us consider system / and forget, for the 
time being, the coupling conditions. We simplify the notation by writing Q instead 



of Q{1). Bearing in mind the coordinatization (3.24), we notice that all the functions 
ciiU) can be represented as 

fx{Q) = Tt{QX) =< QX > 
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by means of a suitable upper triangular matrix X. We introduce two Poisson brackets 
as follows 

{/x,/y}i(Q) = <Q[X,Y]> (3.39) 
{/x,/y}2(Q) = <{XQ)+YQ>-<{QX)+QY> (3.40) 
+ < QVI+Y > - < QYI+V > 

where 

oo j 
j=0 1=0 

In terms of the coordinates we have 

{am{i),an{j)}l = am+n{j)Si,j-n - am+n{i)Sj,i-m 

and 

{am{'i),an{j)}2 = an+m+l{j)Si,j-n-l — an+m+l{i)Si,j+m+l 

m 

n+m—l ij)^i,j—n+l CLn+m—l (0*^^ {j)^i,j+m—l ) 

j 

+am{-i)an{j) ^ ((^ii - Si-m,l) 
l=j-n 

Introducing the Hamiltonians = -f^fc(l) = ^Tr((5'^), one can verify that the two 
Poisson brackets are compatible, i.e. 

{Hk+i,fx}i{Q) = {Hk,fx}2{Q), 

that they are in involution and that they generate the flows with respect to ti i^.. 

With minor modifications we could do the same for system II. For example, in 
the 2-matrix model we can find another couple of compatible Poisson brackets and 
another set of Hamiltonians Hk{2) in involutions. As we already noticed the flows 
generated by the two series of Hamiltonians commute. However we are not able to 
decide whether the -fffc(l)'s commute with the -fffc(2)'s. If not, this would be an 
example of integrable systems with non-commuting hamiltonians. 



4 VFi+cx) Constraints in Two— Matrix Models 



It was shown in |1£] that the path integral of multi-matrix models has a large 
symmetry that can be translated into algebraic constraints on the partition function. 
It is instructive to see how this large symmetry is inherited by the corresponding 
linear systems. We will limit ourselves here to two-matrix models. Starting from the 
string equations and the discrete hierarchies we found in the previous section, we will 
prove that the partition function of the two-matrix model satisfies Wi+oo~constraints. 
To this end it is actually convenient to consider a much more general model than the 
one studied so far, which is characterized by the most general polynomial coupling 
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between the two matrices. The reason is that this enlarged model possesses a larger 
symmetry than the original one. By using the ensuing algebraic structure our goal 
will be reached in a much easier way. Eventually, by reduction, we will find the result 
appropriate to the original 2-matrix model. 



4.1 The generalized coupling conditions 

Let us introduce a more general polynomial interaction term into the partition func- 
tion 

„ 2 Af 

ZN{t,C) = consty n n^V*^(^i)^(^2)e", (4.1) 

a=l i=l 

oo N oo N N 

^ = E*i,'«E^t + E*vE^2.+ E ^'^^En.^. (4-2) 

k=l 1=1 r=l i=l a,b>l i=l 

In order to recover our original two-matrix model we have to simply set Cab = 
except for Cn = ci2, where the latter is the coupling of section 2 |. 

We follow the same procedure we followed in section 2 after integrating out the 
angular part, and introduce the orthogonal polynomials 

e„(Ai,t) =A^ + ---; 7]n{X2,t) = X^ + ---. 

which satisfy the improved orthogonal relation 

d\id\2in{Xi,t)e^ rim{\2,t) = hn5nm- (4.3) 

where 

UiXi,\2) = Vl{Xl)+V2iX2)+ E CabX^l>^i 

a,b>l 

oo oo 
ViiXi) = ti,fcA^ V2{X2) = E t2,rX2- 

k=l r=l 

Next we define the Q matrices 

dXidX2in{Xi,t)e^ Xa'qni{X2,t) = Q„m{a)hm = Qmn{a)hn, « = 1, 2 (4.4) 



and their conjugates 



We appended the label o in order to distinguish these operators from the previously 
defined P(l) and P(2). 

As in section 3 we can get the spectral equations 

Aie = Q(l)e, A27? = Q(2)7?, (4.5) 



§At times, eqs.(4.19) and (4.25), we will use a compact notation in which Co,b = t2,6 and Cafi — ti.a 



14 



and the equations of motion of the polynomials, in the form of two coupled linear 
systems, which coincide with the systems I and //of the previous section when q = 2. 
The only differences can be found when the coupling conditions are involved. Noting 
the spectral equations, one can easily see that the coupling conditions become^ 

Po(l) + Fi(l) + E aCabQ''-\l)QH'^) = 0, (4.6a) 

a,b>l 

Vo{2) + ^2(2) + bCabQ''{l)Q'-\2) = 0. (4.6b) 

a,b>l 

The time evolution of the partition function is given by eq.( 3.30| ), with a = 1 and 
2. We can also write down the coupling dependence of the partition function 

lnZ^(t,C7) = Tr(Q'*(l)Q^(2)), Va,6>l. (4.7) 



dCab 

This is all we need for the derivation of the constraints. 



4.2 Virasoro constraints 

The VFi+cxD-algebraic constraints contain a particularly important subset, the Vira- 
soro constraints. We begin with them. Since the discussion is valid simultaneously for 
both linear systems, we temporarily omit the system indices and consider a general 
Jacobi matrix Q and its conjugate P 

Qij = 6j^i+i + Sj6ij + Lj6j^i-i + ..., [Q,P]=1. (4.8) 

As was shown in [20| the string equation completely determines the matrix P. Pro- 
ceeding as in |2C], we can express the string equation as a series of equations 



Tr(Q"+i(l)(Po(l) + l^'(l)+)+ E aCabQ''~\l)Q\2))]=0, n > -1 (4.9) 

Now using the relation between partition function and Q matrix, these are equivalent 
to 

(£W(l)+/li](l))Z;v(i;C) = 0, n> -1. (4.10) 

where 



oo 



^Hereafter we will use the notations 

oo 

n'(l) = E^*l^^-Q'~'(l)' 

oo 

n"(i) = E^(^-i)^i.'^-^'"'W' 



k=2 

etc. 
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cm 



oo 3 \ 



oo 
k=l 



d n + 1. d 1 



n-l 



k+n 



+ (iV + 



52 



2 5ti,„ 2 dti^kdh^n-k 



and 



rW(i) 



51 0,CabT^ + 



a>2,b>l 



J2 "^'^''^r — 

a,6>l OrOa+n.b 



dCa-l,b ^ 9*2,6 



n > 0, 



These operators satisfy the Virasoro algebras 

= (n-m)/:W^(l), n,m > -1, 

[£W(1),^J.'1(1)]=0, n,m>-l, 
[4i](l),rW(l)] = (n-m)r[^|„(l), n,m > -1, 



(4.11) 
, n > 1. 

(4.12a) 
(4.12b) 



(4.13a) 
(4.13b) 
(4.13c) 



In order to recover the result we are interested in for our original two-matrix model 
we should simply set, in eqs.( 4.10| , 4.12a , 4.12b| ), Cab = except for Cn = C12, where 
the latter is the coupling of section 2. 



4.3 Higher rank constraints 

In the previous section we obtained the Virasoro constraints, which may be referred to 
as rank 2 tensorial constraints. In the following we deal with higher rank constraints. 
In order to get the spin-3 operators, we introduce the following notations 

V = ^ C,,Q"(1)Q^(2), 

a,b>l 

V' ^ aCabQ''-\l)Q\2), 

a,b>l 
a,b,a ,b' >1 

From the trivial relation 

J (iAi(iA2^(en(Ai,t)e^r?„(A2,t))=0, 

we get the following identity 

- + + y'\i^ + 2y'(l)V'+ : V'' : +V" = 0, (4.14) 

We claim that multiplying ( [4.14| ) by Q"'~^'^{1) from the left and taking the trace, we 
obtain the rank-3 constraints 

(£[2](i) _ r[2](l))Z^(t; C) = 0, n> -2. (4.15) 
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with the definitions 

OO o CXD 1 — 3 r)2 

+2iV(E iti,iJ^ + Nh,2 + Itl,), (4.16a) 

OO o OO 1—2 o2 

£5(1) = E hi2h,i,h,i,^ + E^^ME^r^^ 

+2^E^*M^7^ + ^-i(l) + ^V> (4-16b) 

OO rj OO i — 1 r\2 

£!fi(i) = E ^i^2iMi*i-^.^^ — + E^*i>'E 



d 1 

+2(7V + 1) E + ^iV(iV + 1)(7V + 2), (4.16c) 

1=1 ^hi 3 

OO o OO l+n~l o2 

£[?i(i) - E ^i^2ti,,ti,.^-^^ — + E^*M E 



1 n— 2n— 1 o3 

+577^1: EC" -' + 2)- 



2(n + 3) ^ dti^idti^kdti,n-i-k 
,2 . . . .... . 1. - . ... d 



and 



+(Ar^ + (n + l)Ar + -(n + l)(n + 2))- 
+{2N + n + 2)C\l\l) n>l. (4.16d) 



r3(i) ^ E ^^^'CabC,, ^ 



a+a'>3,b,b'>l dC^+a'-2,b+b' 

+ E a{a-l)Cabj-^ 

a>3,b>l 0Ca-2,b 

+ E^-^i.'^ + 2EC2.^, (4.17a) 



^^a+a'+n,6+6' 



a,b,a ,6 >1 

+ E «(« - l)Ca6 ^^^ , n > -1. (4.17b) 

They satisfy the following algebra 

[/:P1(1),£W(1)] = (n - 2m)4%„(l) + m(m + 1)£L'U(1) (4.18a) 

(1), (1)] = 2(n - m)£!fU(l) - (n - m)(n + m + 3)£PU(1) (4.18b) 

[rPl(l),TW(l)] = (n - 2m)rP|„(l) - m(m + 1)4^1^(1) (4.18c) 

[rPl(l),rPl(l)] = 2(n - m)r[=!l^(l) + (n - m)(n + m + 3)ri5!„(l) (4.18d) 
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where we have introduced the spin-4 operator 



(1) 



E 



d 



ai,bi,a2,b2,a3,b3>l 

3 

+ 



dC. 



ai+a2+a3+n, 61 +62+63 



E 



ai,fei,a2,fe2>l 

+ a{a-l){a-2)Cab 



aia2{ai + a2 — '2)CaibiCa2b2-K^ 

^ 0'^ai+a2+n,bi+b2 



d 



a,b>l 



dC, 



n > -3. 



a+n,b 



(4.19) 



As for £n (1), it is a very comphcated expression, which will not be written down 
here. One can use eg. (4. 18b) as its definition. 



Now let us prove our claim — the constraints (4.15). At first, we notice that 

Tv(p2(i))=o, 

So taking the trace of eg. ( [4. 14 ), we can easily rewrite it as follows 

(£L'kl)-r3(l))^Jv(t;C7) = 0, (4.20) 

which is a particular case of egs.( |4.15| ), when n = —2. Then, the algebras ( [4.18aD 
and ( [4.18c ) guarantee that all the other constraints in egs. (|4.15D are true as well. On 
the other hand from ( [4.18b| ) and ( [4.18d ), together with ( [4.15| ), we obtain the rank 4 
operator constraints, 



n > 



(4.21) 



We see that a larger algebraic structure helps us simplifying the calculations dras- 
tically. Therefore it is evident that we should keep introducing higher and higher 
tensor operators in order to close the algebra generated by {(£^^'(1) + Tn\l)),n > 
-1; {C^n\l)-Tl^\l)),n> -2}. For example, from egs.( |418bD and ( |4l8dD , we obtain 
the rank-4 operators £n (1) and Tn (1). In turn the commutators of these operators 
with Cn (1) and TA (1) generate the rank-5 operators and so on. In this way we 
must introduce the whole sequence of operators 

^'''(1) = E aia2 . . . arCa-,bi ■ ■ ■ Ca,br 



d 



dC, 



a\+...+ar+n,b\+...+br 



aia2 

ai,bi,...,ar,br>l 

+ lower orders 

and the corresponding Cn\l)- All these generators form a closed W^i+oo-algebra 



[4'-l(l),rH(l)] = {sn - rm)Tl'Z-'\l) + ... ^ 
[£H(1),£W(1)] = (sn - rm)£|:+r^l(l) + . . . , 
[7l'^l(l),/:W(l)]=0, 



(4.22a) 
(4.22b) 
(4.22c) 



for r, s > 1; n > — r, m > —s. Here dots denote lower than r + s — 1 rank operators. 
Since we have already shown that a combination of the rank-2 and rank-3 operators 
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annihilate the partition function, eqs.( 4.10| ) and ( f4.15 ), then, due to the above H^i+oo 
algebra [J], the same combination of higher rank generators will annihilate the partition 
function, that is to say 

(£W(l)-(-ir4^1(l))^7v(t;C)=0, r>l; n > -r. (4.23) 

In the same way from the second linear system we can derive another copy of 
Wi+oo algebra as well as another version of the VFi+oo constraints 

(/:H(2) - {-lYT}{\2))Zr,{t;C) = 0, r > 1. (4.24) 

where c\^\2) can be obtained from Cn\l) by simply replacing the ti^^ parameters 
with the t2,r ones. However, as for the Tn\2) generators, we have 

4'](2) = J2 bCabjr^. n>-l 

a,b>l Cl'-'a,b+n 

T^\2) = bi;C^,C,, ^ 



'^^ a+a ,h+b' +n 



+ Y^b{b- l)Cab ^^^ , n > -2, (4.25a) 

a,b>l OUa,b+n 

= hh2b-iCa^biCa2b2Ca-ib3j^ 

ai ,a2 ,62 ,a3 ,63 > 1 C'-^ai +a2 +a3 +62 +63+" 

3 d 

+ Y 0^1^2(61 +b2- 2)Ca,biCa2b2- 



ai,fei,a2,fc2>l 



OLyai+a2,bi+b2+n 



-YKb- m - 2)Cabj^p^, n > -3, (4.25b) 

a,b>l OCa,b+n 



rM(2) = Y hb2-..brCa,b,-..Ca^b.. ^ 



aiM,-,ar,br>l dCai+...+arM + -+br+n 

+lower order terms, (4.25c) 
with the following algebra, which is isomorphic to (4.22a- 4.22(^ ) 

[rM(2),rW(2)] = (sn-rm)7;^;+^-^](2) + ..., (4.26a) 

[/:M(2),/:H(2)] = {sn - rm)/:![+r"(2) + • • • , (4.26b) 
[Tl^\2),d^{2)]=0, (4.26c) 

for r,s > 1; n > — r, m > —s. 

In conclusion, we see that, in the case of a general interaction, there are two 
isomorphic VFi+oo algebra constraints . The closed algebra formed by these two 
pieces of VFi+oo algebras gives the complete constraints (one should notice that they 
do not form a direct product). 



"in fact we have two Wi+oo algebras, one formed by the T generators and the other by the C generators, 
but since for the constraints the relevant algebra is the direct sum r„ + Cm we keep speaking about a 
unique algebra. 
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4.4 VFi+oo constraints 



From eq.( |4.23 ), it is easy to see that if we return to our original model, i.e. we set 



all the Cab, a,b > I equal to zero, and only keep Cn = C12 = c 7^ 0, we have 

TM(l)Z^(t;C) = c'^Tr(Q^+"(l)Q'-(2))z^(t;C) (4.27a) 

TM(2)Z^(t;C7) = c^Tr(Q'^(l)Q'-+"(2)) Z;v(t; C). (4.27b) 
Substituting these into eqs. (|]2|) and ( PI , we get the following constraints 

wWZjv(t, C) = 0, r > 0; n > -r, (4.28) 

where 

T^M ^ (-c)"£M(l) - (4.29) 

form a algebra. One can explicitly check that this result coincides with the 

one in ref.[19|. 

A further reduction is possible. Suppose we set t2,fc = 0,k > q, then from the 
above equation, we have 

&i{l) + d'TT{Q''{2))=0 (4.30) 

Substituting it into the other constraints, we can get other VFi+00 constraints, which 
are only expressed in terms of t2,fe with 1 < k < q (beside all the ti,fc's, of course). 
The operators form a subalgebra of the Wi+00 in eq.( 4.2^ - ^.29| ). 



5 Differential Hierarchies of Multi— Matrix 
Models. 

So far we have shown that multi-matrix models can be represented by means of 
coupled discrete linear systems, whose consistency conditions give rise to discrete KP 
hierarchies and string equations. In this section, we will use the method introduced in 
[^] to transform the discrete linear systems into equivalent differential systems whose 
consistency conditions are purely differential hierarchies. In the most general case 
they coincide with the generalized KP hierarchy. In this elaboration no continuum 
limit is involved. 

As anticipated in the introduction the clue to the construction are the ti^i and 
tq^i flows. On the one hand, it is just the first flow equations of ^ and ^ that enable 
us to eliminate the difference operations from the RHS of the other flow equations 
and to recast the discrete linear systems into a purely differential form. On the other 
hand, using them, we can express the flows of the coordinates of the j~th sector as 
functions of the coordinates in the same sector. This remarkable property is far from 
obvious a priori. 

{*) Differential linear system I. 

We apply the procedure outlined above to the discrete linear systems derived 
in the previous section. Since for q-matrix models we derived 2q discrete linear 
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systems, when we pass to the differential language we expect to find 2q differential 
linear systems. We start with the DLS (| 
have 



Using the eqs .( 13.341 ) and ( 3.33 ), we 



ao(j - 1) = ao{j) - {InRjY 
aoU - 2) = ao{j) - {InRj)' 

where for any function f{t), we write 

df 



In 



ao(j) + ^ 



du 



In general, defining 



R 



■J+r 



Of, 



Rj^r=F-{j), 



df 



«o(j - r) 



df 



G-{j) 



we obtain the recursion relation 



F++i(i)=F+(j)+G+(i) 

G++i(j) = G+(i) + (ln[F+(i) + Gtm 



F;_^{j) = F-{j)-G-{j) + 

G-_i(j) = G-(j)-(lnF-)' 



dti idt, 



9,1 



(5.1a) 
(5.1b) 

(5.1c) 
(5.1d) 



These results guarantee: 1) that all the ao(i)'s and i?j's {i ^ j) can be expressed in 
terms of ao{j) and Rj and their ti^i and tg^i-derivatives; 2) that, substituting these 

results into eqs.( |3.27a - 3.27f ), we can recursively express all a;(i)'s and Tl^"\iys{i ^ j) 
as functions of the coordinates in the j-th sector. 

Now let us consider the ti^i-,and tg^i-flows of ^j. From eqs.(^), we can write 
down their explicit forms 



d 



-Ri^ 



which lead to the following equalities 



(5.2) 
(5.3) 

(5.4) 



where 



Aj = d- ao{j) = -d ^Rj+i 

1=0 



R 



-d 
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It is easy to see that 



1 Vj > 1 



Using eq.(^.4|) we can rewrite the spectral equation in 
equation 



where 



Lj{l) = d+ ai{j)Bj^iBj^i+i . . . 
1=1 



purely differential 
(5.5) 



1 



1 



CO 



d-ao{j-l) d-ao{j-l + l) d - ao{j - 1) 



1=1 



(5.6) 



where [j, mi] denotes the least between j and mi. All the functions ui{j) are functions 
of the coordinates in j-th sector only. 



U2U) = a2{j) + ai{j) ^ao(j) - (InRj)' 
U3{j) = 03 (i) +a2{j) 



(5.7) 



2ao(i) - 2{lnRjy - { In R, - ho{j) + 



+ ai(i) [aoU) - {InR,)']' - a',{j) + (Inii,)" 



and so on. Lj{l) is an operator of KP type. Actually it is in general a reduction of 
the KP operator 

00 

LKP = d + J2wiij)d-' 

1=1 

where wi are generic (i.e. unrestricted) coordinates. We notice that due to the 
limitation in the first two summations in ( |5.6D , the operators Lj(l) do not have in 
general a universal form in terms of the coordinates ai{j). In fact they change with 
j as long as j < mi, while the form of Lj{l) is independent of j if j > mi. We will 
comment upon this point later on. 

Without repeating all the steps of |13], the rules to pass from lattice to purely 
differential language are: 

(i). The Jacobi matrices are mapped into KP-type operators, i.e. 



(^Q{a)^^,^Lj{a)^j, a = l,2,...,g 
(P(1)^).^M,(1)^',. 



(5.8) 
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where 



L(a) 



1=1 



(5.9a) 



^'o^^ (i) + E (i)4+/-i^+^-2 • • • i.- 
+ E^/'^(i)4-/4--w-4-i 



(5.9b) 



We remark again the non-universaUty of such expression except for j large enough. 

{ii). The lower triangular part of a Jacobi matrix (or its powers) maps to the 
pure integration part of the KP operator (or its powers); 

(Hi). The upper triangular part together with the main diagonal line of the Jacobi 
matrix (or its powers) correspond to the purely differential part of the KP operator 
(or its powers). 

(iv). The residue of the KP operators have a particularly simple form, 



resd{Lj{l)) = ai{j) = Qjj-i(l). 
More generally, we can obtain 

and 



(5.10) 
(5.11) 



as well as 



ves9iLj{a))=Tt\j)=Qj,j-iia). 



res 9{ma)) = {Q'{a)y,,Ja) 



(5.12) 



Expanding all the above operators in powers of d, the coefficients can be expressed 
as functions of the coordinates in the j-th sector. Collecting all the results we get a 
linear system 



-(L5(a))_* 



9 ^ 

Otlr J 



a = 2,3,...,q 



(5.13) 



The consistency conditions are 
d 



dtpr 

d 



dt 



L,{a) = [mm^,Lj{a)], 
L^{a) = [Lj{a),mm_], 



l<P<a 
a<p<n 



(5.14) 
(5.15) 



23 



We do not drop here the label j in the above equations in order to stress the fact 
that, although the hierarchy is expressed in terms of the coordinates of the j-th 
sector only, Lj depends on j, as we have already remarked. 

The linear system (|5.13| ) is exactly the same as the one obtained in [|l5| in a 
completely different context, i.e. by trying to generalize the KP hierarchy. 

We remark that if we impose the condition 

aiU) =0, V / > 2 

then the second expression of (|5.6| ) gives the so-called two-bosonic representation of 



the KP hierarchy. Motivated by this fact, we refer to the full expression of eq.(5.6) 
as the multi-bosonic representation of the KP hierarchy. 

(**) Differential linear system II. 

In the above analysis, we considered as the space coordinate, and d as the 
basic derivative. However, as we already pointed out, it is also possible to consider 
tq^i as the space coordinate and d as the basic derivative. In this way we would get 
another linear system, in which <I> plays the role of Baker-Akhiezer function. In the 
following we discuss this system. As before we have 

^.+1 = C,'^,, $,-1 = b,-i^, (5.16) 

where 

Cj=d-ho{j) = -d-^Rj+i 

oo -. 



1=0 ^^+1 



The linear system is 



Lj(5)<I>j = A,$j (5.17a) 

= -(Z;Xa))_$„ a = l,2,...,g-l (5.17c) 
Mj'^j = (5-17d) 

Now the KP operator takes the following form 

j 

Ljiq) =d + J2 bi{j)bj-ib,^i+^ . . . b,.i 
1=1 

= d+ybi{j)^ — - — — - — 

oo 

= d + J2viij)d-' (5.18) 
1=1 



24 



All the functions vi{j) are only functions of the coordinates in the j-th sector. 

vi{j) = bi{j), 



V2U) = b2{j) + hU) boij)-d{lnRj) 



(5.19) 



V3U) = Hj)+b2U) 



26o(i) - 2d{lnR,) -dllnR,- dbo{j) + 



+ bi{j) Mj) - d{lnRj)r - dboij) + d'lnR, 



d'^lnRi 



dti idt, 



9,1 



and so on. Using the rules to go from the lattice to the differential language, we 
obtain 



(Q(a)$) 
(P(,)$)^. 

with new differential operators 



on 



a = 1,2, . . . ,g 



(5.20a) 
(5.20b) 



1=1 

00 

MM) = H'Hj) + E Pi'hj)Cj+i-iC,+i.2 



(5.21a) 



+$:p/'')(j)i3,_,^,_,+i.. 

1=1 

Therefore we can write down another linear system 



D 



(5.21b) 



$0 



d 



(5.22) 



a 



1,2,...,, 



I M,(l)$, 



whose consistency conditions also give the hierarchy ( |5.14 " 5TT5| ). 

We could also start from the general form of the discrete linear systems ( ^.22 ) 
and ( |3.23| ), translate them into the differential language by using the rules we listed 
before. Once again they lead to the hierarchy 05.14 - 51151 ). We will spare the reader 
the details. 

We end this section with a remark on the flows in ( 5.14 ) and ( ^.15| ). They all 
commute with one another, since their lattice versions do. So all the differentiable 
hierarchies ( |5.14| ) and ( ^.15 ) are integrable. This is a very remarkable result. We 
know in fact that attempts at enlarging the KP hierarchy in a very general way 
have led to incompatible KP flows 1 24]. Therefore we may say that the multi-matrix 
models give a natural realization of commutating KP flows. 
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5.1 The coupling conditions 

Our final task is to reexpress the coupling conditions ( |3.5a - 3l5b| ) in the differential 



language. This can be done very easily by replacing matrices by their operator ver- 
sions. The only problem we should care about is the matrix V{q), since its differential 
version is not M(q). We may simply denote its differential version by M{q), then 
the coupling conditions are 

M(l) +ci,2L(2) = 0, (5.23a) 
Ca-i,aL{a - 1) + V'{a) + Ca,a+iL{a + 1) = 0, 2<a<q-l; (5.23b) 
Cq-i,qL{q-l)+M{q)=0. (5.23c) 

where for simplicity we dropped the label j, denoting the sector. Of course 

r=l 

Some of these equations are also met in |15], so they are characteristic of the gen- 



eralized KP hierarchy. But not all of them appear simultaneously in the generalized 
KP hierarchy. More precisely, if we are in the ti^i-picture, that is ti^i is assumed to 
be the 'space' variable, then in the generalized KP hierarchy there naturally appear 
only ( ^.23a| ) and ( 5.23b| ); therefore in this pictu re (|5.23c ) is a constraint on the gen- 



eralized KP hierarchy. Vice, in the t^^i-picture ( 5.23c| ) and ( ^.23b| ) appear naturally 



in the generalized KP hierarchy, while ( ^.23a| ) plays the role of constraint. 

Thus we may conclude that multi-matrix models correspond to the generalized 
KP hierarchy subject to certain constraints. 

5.2 The partition function as r— function 

In the generalized KP hierarchy we can define the following r-function 



dUidtc 



-Inr = les qL^ (a), V a,r (5-24) 



where L is the generalized KP operator. We want to prove that the partition function 
of multi-matrix models is a r-function of the a-th generalized KP hierarchy. To 
achieve this we only need to prove that eqs.( |3.3i| ), once translated into the differential 
language, have the same form as eqs.( 5.2^ ). This can be easily done, if we note that 



eq.( 5.12 ) is valid for any 1 < a < q, and any positive integers r, j. Choosing in 



particular j = N, we have 

This equality together with eq.( |3.3l|) gives 



dtiAdto 



-lnZ]y{t,c) = res g(^L]\f{a)j (5.26) 



This is nothing but ( ^.24 ). Therefore we can conclude that the partition functions 



of multi-matrix models are exactly r-functions of the generalized KP hierarchy. 
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6 Examples and Conclusion 



In the previous sections we were able to extract differential hierarchies from multi- 
matrix models, without passing through a continuum limit but by simply studying 
the properties of the matrix model lattice. These arc KP type hierarchies of which 
the original partition function is a r-function. To determine the initial conditions 
for this T-function we digged out additional equations: the coupling conditions or 
string equations. 

What we have worked out so far is a very general scheme which allows for a vast 
number of possibilities, i.e. physical models, the physical meaning being attached to 
the coordinates we have introduced. The analysis of all the different possibilities is 
out of the question. Our present purpose is much more modest: we want to single 
out significant subclasses of models. To do this we can play with two tools: the first 
is the type of potentials introduced in the original multi-matrix path-integral; the 
second tool, once the potentials are fixed and the corresponding differential linear 
systems are worked out, is to study the various possible reductions of the latter. 

To start with we want to split all the possibilities wc arc faced with in two broad 
classes or alternatives. Wc recall what we repeatedly pointed out in the last section: 
the differential hierarchies obtained from the lattice are in general not universal in 
that they depend on j, i.e. on the sector. In particular, the flows of the logarithm of 
the partition function depend in general on the sector. For the differential hierarchies 
coming from a given lattice hierarchy we have therefore two possibilities, which are 
connected with specific choices of the numbers Pa characterizing the order of the 
potential V^: 

- the non-universal alternative: the differential hierarchies coming from the same 
lattice hierarchy are all different in different sectors. In this case one can hardly 
consider these hierarchies an intrinsic property of the lattice. Presumably the only 
way (if any) to treat this case is via a continuum limit. 

- the universal alternative: the differential hierarchies coming from the same lat- 
tice hierarchy become isomorphic, i.e. truly independent of j, for j large enough. No 
doubt, the j-independent hierarchies, referred to henceforth as the relevant univer- 
sal hierarchy, express intrinsic properties of the lattice, and they can presumably be 
realized in terms of topological field theories. Needless to say, we are interested in 
this case. 

Let us see, as an example, the two-matrix model case. We have two parameters: 
Pi and p2- On the basis of the discussion in subsection 3.1, we see that the lattice 
system I gives rise to a universal hierarchy only if p2 is finite, while the lattice system 
//does the same only if pi is finite. Generalizing this discussion to other multi-matrix 
models is straightforward. 

This distinction between universal and non-universal alternative is perhaps diffi- 
cult to grasp at first reading. However it has a rather natural meaning: we recall once 
again that the information concerning the multi-matrix model is contained in the 
partition function; the fiows of the latter depend on the sector j = N (see subsection 
5.2); in the universal alternative the potentials are simple enough as to allow us to 
retrieve all the information contained in the lattice if we choose a sufficiently large 
(but finite) N, while of course only partial information is retrieved if N is too small; 
in the non-universal case the potentials are too complicated for us to retrieve all the 
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lattice information with a finite N, but we might succeed perhaps when N oo.** 
To end this paper we would like to give a few more explicit examples of the 
differential systems of the last section. We will consider the differential system I 
in multi-matrix models; not in full generality however, only in the dispersionless 
This limit is in any case very interesting. It represents the genus 



limit |§, |2f 

approximation and it might characterize the full hierarchy. So let us concentrate on 
the system I of the previous section and, in particular, on the differential operator 
Lj{l), eq.(^), with the following specifications: pi = oo and p2 finite, so that Lj{l) 
satisfies a universal hierarchy for large enough j. For this reason we will drop in the 
following the label j. 

The technique to obtain the dispersionless hierarchy is very simple once we know 
the full hierarchy: in the flow equations we simply discard terms containing more 
than one first order derivative or higher order derivatives. A simple way to obtain 
it without previous knowledge is to replace d with a variable p which is Poisson 
conjugate to ti, 

{P,ti} = 1 

and fully exploit the integrable structure. Therefore we start with the replacement 

1 



1=1 



(6.1) 



where 



Ul+l 



E 

k=0 



^ ) ak+iS'' ^ 



ai 



for / > V 



In ( |6.1j ) we set oq = S in analogy with ||T^. At this point, in order to find the fiows 
of these coordinates, we could simply use the consistency conditions of the system I 
where commutators are replaced by corresponding Poisson brackets. A better way is 
to use the Poisson algebra. The first Poisson bracket is given by 

{a^{x),S{y)}i=5'{x-y) (6.2) 
{ai{x),aj{y)}i= (i + j - 2)ai+j_25 + (j - l)a-+j„2 '^(x - y), 2<i,j<u 

while all the other brackets vanish. The second Poisson bracket is 



{Six),Siy)}2 
{S{x),ai{y)}2 
{Six),aj{y)}2 

{ai{x),ai{y)}2 



!^6'{x-y) 
dS6{x - y) 
!^da,.,5ix - y) 

aid + dai 6{x — y) 



(6.3) 



**We advise the reader not to rely too much in this regard on the analogy with 1-matrix models. The 
analogy between 1-matrix models and multi-matrix models is only partial. One-matrix models are 
all to be classified in the universal alternative, as their Jacobi matrix has at most three non-vanishing 
diagonal lines. The dependence of equations and differential operators on the sector is instead present in 
1-matrix models as well, although in a quite unconspicuous way as it appears only in the zeroth sector 
while everything takes the universal form starting from the first sector on. 
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{ai{x),aj{y)}2 
{ai{x),aj{y)}2 



(i + l)ajd + ja'j'^dix - y) 
(^iai+j_i5 + jdai+j-2 + {i - l)ai+j^2dS 

i-2 

+E 

1=1 



i+j-2 

I - l)ai+j_i_2dai + {j - I- l)aidai+j_i^2 



In these equations i,j > 2. Next we write down the Hamiltonians 

TLr 



(-1) 



5{x - y) 



(6.4) 



The RHS denotes the coefficient of in . The two Poisson brackets are compat- 
ible with respect to these Hamiltonians. Using this fact we can easily write down the 
flows in terms of 



5S ' 



Gr- 



5ai 



5a2 



and recursion relations for the latter. For example, in the case of three flelds, i.e. 
ly = 2, we obtain 

dS _ 



dai 
dtr 



K+i 



(6.5) 



— Za2rjr+1 + 0'2^r+l 



The recursion relations are 



(6.6) 



G;+i = -F; + (5G,)'+-(aii?,)' 

= SFl. + 2aiG', + a[Gr + 3a2E[. + 2a'2Er 
2a2-B;,+i + a'2Er+i = ^aiF/ + SaaC; + a'^Gr 

+ (2a25 - laj)E'^. + (asS 

2 1 

the initial conditions being Fi = Ei = Gi = 0. 

Prom this result we can learn something important. In |13],|14| we obtained the 
KdV hierarchy reducing the NLS hierarchy. The latter is generated as the consistency 
condition of system I for a KP operator of the type L(l) with only two non-vanishing 
coordinates uq = R and S. The reduction to the KdV hierarchy is obtained by setting 
5 = 0. Therefore it would seem not unmotivated to expect that setting 5 = in (|6.5|) 
one finds the dispersionless 3-KdV (or Boussinesq) hierarchy. However an explicit 
calculation proves that this is not the case. In other words this approach to find 
the n~KdV hierarchies as reductions of the differential systems ensuing from multi- 
matrix models is too naive. A way to obtain the n-KdV hierarchy could be as follows: 
reduce L{1) in just the same way we reduce the KP operator down to the n-KdV 
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operator; if the order of the V2 potential (in two-matrix model) is large enough, we 
can obtain all the flows we wish of the n-KdV hierarchy. 

In any case it is evident that the reduction to the n-KdV hierarchies is a very 
drastic one: most of the generalized KP structure gets lost in this passage. Taught 
by the 1-matrix model example, |14], we think a lot of information, in other words 
of topological models, is contained in the generalized KP structure as it appears in 
multi-matrix models. We intend to deal elsewhere with this part of the analysis of 
multi-matrix models. 
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